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We analyze the possibility of measuring the single spin asymmetry (SSA) in open charm pro- 
duction by linearly polarized photons in the planned E160/E161 experiments at SLAC using the 
inclusive spectra of secondary (decay) leptons. In leading order pQCD, the SSA in the azimuthal 
distribution of the charged decay lepton is predicted to be about 0.2 for SLAC kinematics. Our 
calculations show that the SSA in the decay lepton distribution is well defined in pQCD: it is 
stable both perturbatively and parametrically, and practically insensitive to theoretical uncertain- 
ties in the charm semileptonic decays. Nonperturbative contributions to the leptonic SSA due to 
the gluon transverse motion in the target, heavy quark fragmentation and Fermi motion of the 
c-quark inside the D- meson are predicted to be about 10%. We conclude that measurements of 
the azimuthal asymmetry in the secondary lepton distribution would provide a good test of the 
applicability of pQCD to open charm production at energies of fixed target experiments. Our 
analysis of the SLAC experimental conditions shows that this SSA can be measured in E160/E161 
with an accuracy of about ten percent. 
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I. INTRODUCTION 

In the framework of perturbative QCD, the basic spin-averaged characteristics of heavy flavor hadro-, photo- and 
electroproduction are known exactly up to the next-to-leading order (NLO). During the last ten years, these NLO 
results have been widely used for a phenomenological description of available data (for a review see §). At the 
same time, the key question remains open: How to test the applicability of QCD at fixed order to heavy quark 
production? The problem is twofold. On the one hand, the NLO corrections are large; they increase the leading order 
(LO) predictions for both charm and bottom production cross sections by approximately a factor of two. For this 
reason, one could expect that higher-order corrections, as well as nonperturbative contributions, can be essential, 
especially for the c-quark case. On the other hand, it is very difficult to compare pQCD predictions for spin-averaged 
cross sections with experimental data directly, without additional assumptions, because of a high sensivity of the 
theoretical calculations to standard uncertainties in the input QCD parameters. The total uncertainties associated 
with the unknown values of the heavy quark mass, mg, the factorization and renormalization scales, fj-p and fiR, 
Aqcd and the parton distribution functions are so large that one can only estimate the order of magnitude of the 
pQCD predictions for total cross sections at fixed target energies 1^,^. 

In recent years, the role of higher-order corrections has been extensively investigated in the framework of the soft 
gluon resummation formalism. For a review see Ref. B. Formally rcsummed cross sections are ill-defined due to the 
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Landau pole contribution, and a few prescriptions have been proposed to avoid the renormalon ambiguities • 
Unfortunately, numerical predictions for the heavy quark production cross sections can depend significantly on the 
choice of resummation prescription |^ . Another open question, also closely related to convergence of the perturbative 
series, is the role of subleading contributions which are not, in principle, under control of the resummation procedure 

i|. 

For this reason, it is of special interest to study those observables that are well-defined in pQCD. An nontrivial 
example of such an observable is proposed in [ p^[p^ ] , where the charm and bottom production by linearly polarized 
photons, 

7UiV^Q + X[Q], (1.1) 

was considered0. In particular, the single spin asymmetry (SSA) parameter, Aq(jiqt), which measures the parallel- 
perpendicular asymmetry in the quark azimuthal distribution, 

—(pqt^^q) - —-^{pQT)[l + AQ{pQT)V^cos2ipQ], (1.2) 



dpQTdtpQ 271 dpQT 

where 



A (p J.) ^ ^ d^^Q(pQT, V?Q ^ 0) - d^agipQT, V?q = tt/2) 
^ ^'^ T't, d2(TQ(pQT, = 0) +d2(TQ(pQT,V3Q = 7r/2)' 



da 



has been calculated. In (|L2D and (|lj), is the unpolarized cross section, d^aqipQT, (fq) = d^|^^(PQT, Vq), 

is the degree of linear polarization of the incident photon beam and ipq is the angle between the beam polarization 
direction and the observed quark transverse momentum, pqt- The following remarkable properties of the SSA, 
Aq(j)qt), have been observed [ p^ : 



The azimuthal asymmetry (1.3) is of leading twist; in a wide kinematical region, it is predicted to be about 



0.2 for both charm and bottom quark production. 

• At energies sufficiently above the production threshold, the LO predictions for Aq{pqt) are insensitive (to 
within few percent) to uncertainties in the QCD input parameters. 

• Nonperturbative corrections to the 6-quark azimuthal asymmetry are negligible. Because of the smallness of 
the c-quark mass, the analogous corrections to Ac{pqt) are larger; they are of the order of 20%. 

In Ref. [O], radiative corrections to the integrated cross section, 



1 unp / r'\ 

:f^iS,VQ) = ^^-^[l + AQiS)V,cos2^Q], (1.4) 



have been investigated in the soft-gluon approximation. (In Eq. (1.4), \fS is the center-of-mass energy of the process 



(1.1)). Calculations indicate a high perturbative stability of the pQCD predictions for Aq{S). In particular. 



• At the next-to-leading logarithmic (NLL) level, the NLO and NNLO predictions for Aq{S) affect the LO 
results by less than 1% and 2%, respectively. 

• Computations of the higher order contributions (up to the 6th order in as) to the NLL accuracy lead only 
to a few percent corrections to the Born result for Aq{S). This implies that large soft-gluon contributions to 
the spin-dependent and unpolarized cross sections cancel each other in Eq. (|l.3|) with a good accuracy. 



^The well-known examples are the shapes of differential cross sections of heavy flavor production which are sufliciently 
stable under radiative corrections. 
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Therefore, contrary to the production cross sections, the single spin asymmetry in heavy flavor photoproduction 
is an observable quantitatively well defined in pQCD: it is stable, both parametrically and perturbatively, and 
insensitive to nonperturbative corrections. Measurements of the SSA in bottom photoproduction would provide 
an ideal test of pQCD. Data on the D-mcson azimuthal distribution would make it possible to clarify the role of 
subleading twist contributions [icIJl^ . 

Concerning the experimental aspects, the azimuthal asymmetry in charm photoproduction can be measured at 
SLAC where a coherent bremsstrahlung beam of linearly polarized photons with energies up to 40 GeV will be 
available soon ||l3|| . In the approved experiments E160 and E161, charm production will be investigated using the 
spectra of decay muons: 

7^ +7V -> c + X[c] ^ /i+ +X. (1.5) 

In this paper, we analyze the possibility to measure the SSA in heavy quark photoproduction using the decay 
lepton spectra. We calculate the SSA in the decay lepton azimuthal distribution: 



dV, , , 1 da""P 
— (Per, (fie) = — j" 



{peT,fe) = — -JT — {pit)[^ + Ai{piT)'PjCos2(pi], (1.6) 



where (pi is the angle between the photon polarization direction and the decay lepton transverse momentum, pix- 
Our main results can be formulated as follows: 

• The SSA transferred from the decaying c-quark to the decay muon is large in the SLAC kinematics; the ratio 
Ae{pT)/Ac{pT) is about 90% for pr >l GeV. 

• pQCD predictions for AiijiiT) are also stable, both perturbatively and parametrically. 

• Nonperturbative corrections to Ai^pex) due to the gluon transverse motion in the target and the c-quark 
fragmentation are small; they are about 10% for pgT >1 GcV. 

• The SSA in Eq. (^^) depends weakly on theoretical uncertainties in the charm semileptonic decaysi, c 
l^veXq {q = d,s). In particular, 

— Contrary to the the production cross sections, the asymmetry Ai{p£t) is practically insensitive to the 
unobserved stange quark mass, rug, for per >1 GeV. 

— The bound state effects due to the Fermi motion of the c-quark inside the _D-meson have only a small 
impact on Ae{peT), in practically the whole region of per- 



We conclude that the SSA in the decay lepton azimuthal distribution (1.6) is also well-defined in the framework 
of perturbation theory and can be used as a good test of pQCD applicability to open charm production. 

The paper is organized as follows. In Section 2 we analyze the properties of the decay lepton azimuthal distri- 
bution at leading order. We also give the physical explanation of the fact that pQCD predictions for the SSA are 
approximately the same at LO and at NLO. The details of our calculations of radiative corrections are too long 
to be presented in this paper; they will be reported separately in a forthcoming publication p5[ |. In Section 3 we 
discuss the nonperturbative contributions to Aglpix) caused by the Peterson fragmentation, the fej^-kick and the 
Fermi motion of the heavy quark inside the heavy hadron. Section 4 contains experimental considerations. We 
discuss the conditions needed to measure the SSA at SLAC, estimate background contributions and expected errors. 



^For a review see Ref. 
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II. PQCD PREDICTIONS FOR SSA 



A. Partonic Cross Sections 



At the Born level, the only partonic subprocess which is responsible for the reaction (1.5) is the heavy quark 
production by photon-gluon fusion, 



-f\k^) + g{kg) -> Q{pq) + Q{pq) ^ l{pt) + vt + q + Q, 



(2.1) 



with subsequent decay c — > t^v^q [q — d, s) in the charm case and b — > l^Viq {q = u, c) in the bottom one. To 
calculate distributions of final particles appearing in a process of production and subsequent decay, it is useful to 
adopt the narrow-width approximation, 
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(2.2) 



with Tq the total width of the heavy quark. Corrections to this approximation are negligibly small in both charm 
and bottom cases since they have a relative size 0{Tq /mq). 

In the case of the linearly polarized photon, the heavy quark produced in the reaction (|]^) is unpolarized. For 
this reason, the single-inclusive cross section for the decay lepton production in (|2.lD is a simple convolution: 



(2.3) 



<i^Pe Tq J Eq d'^pQ d'^pe 

At leading order, 0{aemCts), the (pQ-dependent cross section for heavy flavor production. 



d^PQ 



(pq) 



2sd^aQ 



(s,ti,Ui,(pQ) 



duidtidipQ 

— [Bq (s,ii,ui) -I- ABq {s,ti,ui)VjCos2ipQ] , 



is given by @|6| 



Bq (s, ti, Ui) = TTeQaemas 



Ml h 



5{s + tl+ Ml) 



(2.4) 



(2.5) 



ASq (s, ii, ui) = ire^Qaemas 



tiUi 



6 {s + ti + ui) 



(2.6) 



with eg the quark charge in units of electromagnetic coupling constant. The corresponding partonic kinematical 
variables are 



Ml = {kj — Pq)^ — mq = —Eq^/s{1 — acos^ 
h = {kg -pqf - mq = -Eq^(l + a cost 



— {k^ -j- kg^ , 



a = Jl-ml/El, 



with 6q and Eq the heavy quark polar angle and energy in the 75 center-of-mass system. 
At the tree level, the invariant width of the semileptonic decay Q — s- (viXg can be written as 



Eed^Tsi 
d^pi 



{x) = Isl{x) 



\VckmY 



(27r)4 



l-x 



6(1- 



(2.7) 



(2.8) 
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Here Vckai denotes the corresponding element of the Cabbibo-Kobayashi-Maskawa matrix, Gp is the Fermi con- 
stant, S — mq/mq and 



(1 - acosOiQ), 



(2.9) 



where Ei is the lepton energy and 9iQ is the angle between the lepton and heavy quark momenta in the 75 
center-of-mass system. 



To perform the angular integrations in Eq. (2^), we use the cosine theorem: 

cos 6q = cos 9{ cos 6iQ — sin 9i sin 9iq cos ipiQ , 



(2.10) 



where ipiQ is the azimuth of the lepton momentum in the frame where the third axis is directed along the heavy 
quark momentum. Substituting dj^flq — s-d^rj^g —dcos diQdipeQ in Eq. ( ^.3[ ), we obtain the following master formula 
for the decay lepton azimuthal distribution: 



1 



v/i/2 



nsT, 



Q 



\PQ\<iEQ 



1 2tt 

dcos0eQls\{x) / d(peQ 



(2.11) 



X ^BQ{ti,ui)+ ( 1- 
where the limits of integration are given by 

cos 9iQ = max 

and 



2sin^ 



I sm 



^ 1 L rnlil-S') 
' a [ 2EiEQ 



ABq (<i, ui) Vj cos2ipi 



\J\ — max2[a, 0] ' 



4i;2-m^(l-j2) 
4i;2 + m2(l-52)- 



(2.12) 



(2.13) 



Note that Eqs.( 2.11 )-( 2.13 ) are applicable not only at the Born level, but also in all those cases when radiative 
corrections have a factorizable form (2.3). At the tree level, Eq.(2.11) can be simplified due to ^-function in the 
Born cross sections (2.4) and (2.5). 



B. Hadron Level Results at LO 



At the hadron level, the overall one-particle inclusive kinematical invariants of the reaction (1.5) are defined as 



U — {k^ — Pi) — —EgVzS {1 — cos 6 f ), 
T^{PN-Pif = -^^^^(l-t-cos^, 

The hadron level cross section has the form of a convolution: 

1 



S = {ky +pn) , 
kg = zpN- 



£;,^(5,r,t/,^,) 

d'^pi 



dzg{z,fiF)^^^-^ (zS, zT,U,(pi>) 
d-^pi 



(2.14) 



(2.15) 



in which g{z, ^p) is the giuon distribution function and /ii? is the factorization scale. The lower limit of integration 
in Eq. ( ^.15 ) can be found from the kinematical restriction x < 1 — 5^ . The result is: 



l-max2[/3„,0]' 



(2.16) 
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where 



U 



1 - 



\ 



u 



4UT 



,(1-^2) ) m2 5(1 -52)2 



(2.17) 



Experimentally, to suppress the background contributions originating from semileptonic decays of light hadrons 
and Bcthc-Hcitler process, a cut of the lepton energy in the reaction ( |l.5[) is usually used jl^. For this reason, 
apart from the usual differential distribution. 



dV, 
dperdip. 

we consider also the quantity 

dV, 



iS,T,U,ipe) 



X-dE* 



2 a/ £"^2 _ p2^ 



(2.18) 



Eed^ae 
d^pe 



{S, T, U, (fit) 



\-dE*, 



dpirdipi'' t.cui.x-^^ > r-/ J ^i^^ V. , 2^Ef- ~ pjj. 

In Eqs. ( |2.1^ ) and ( pls| ), E^ is the lepton energy in the lab (nucleon rest) frame, 

with the maximum value 



(2.19) 



2mM 



(^X+VS{1 + /3)(1 - 5^) + A_^5(l + /3)2(1 - <52)2 _ iQp2^^ 



(2.20) 



(2.21) 



and an experimental cutoff E'^ 



,cut' 



A. 



1 ± to2,/5' and P = Jl - im^/S 



Let us discuss the hadron level pQCD predictions for the asymmetry in azimuthal distribution of the decay 
lepton. In this paper, we will discuss only the charm photoproduction at the SLAC energy E^ w 35 GeV, with 



E^ ^ {S~ 



)/2mN- Unless otherwise stated, the CTEQ5M [h7l parametrization of the gluon distribution function 



is used. The default value of the charm quark mass is mc = 1.5 GeV. 

Our calculations of the quantities A^{pt) and Ac{pt) are given in Fig.|^ by solid and dashed lines, respectively. 
One can see that the asymmetry transferred from the decaying c-quark to the decay muon is large in the SLAC 
kinematics; the ratio A^(pt)/^c(pt) is about 90% for pT >1 GeV. Note that pT = pqt when we consider the heavy 
quark production and pt = Pit when the quantity A^(pit) is discussed. 

In Fig. H we show the LO predictions for A^{pt) at different values of the cutoff parameter -E^j-uf One can see 
that, contrary to the production cross sections, the SSA depends weakly on i?|j.ut ^^^^^ practically the whole region 
of Pt- This property of A^{pt) will be especially useful in interpretation of data. In Fig. || we show the SSA, 
■^ij.{pt)i in different intervals of the muon energy E^. 

The most interesting property of the azimuthal asymmetry, closely related to fast perturbative convergence, is its 
parametric stabilityO. As was shown in llo| ], the LO predictions for Aq{pt) are insensitive (to within few percent) 
to standard theoretical uncertainties in the QCD input parameters: ^r, fip, ^qcd and in the gluon distribution 
function. We have verified that the same situation takes place in the case of decay leptons too. In particular, all the 
CTEQ5 versions of the gluon density, as well as the MRST |l^ parametrizations, lead to predictions for Ap(pj') 
that coincide with each other with an accuracy of better than 1-2%. 



^Of course, parametric stability of the fixed-order results does not imply a fast convergence of the corresponding series. 
However, a fast convergent series must be parametricaly stable. In particular, it must be fiR- and /iF-independent. 
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The main source of uncertainties in the LO predictions for both Ac{pt) and Afj_{pT) is the charm quark mass. The 
dependence of pr-spectra on the variation of rUc is shown in Fig. |^ and Fig. || for Ac{pt) and A^(pt), respectively. 
One can see that changes of the charm quark mass in the interval 1.3 < rric < 1.7 GeV affect the quantity Ac{pt) 
by less than 20% a.t px < 2.5 GeV. Analogous changes of rric lead to 20% variations of A^{pt) at 1 < pr < 2.5 
GeV. 

We have also analyzed the dependence of the SSA in the Icpton distribution on the unobserved strange quark 
mass, rris. Fig. ^ shows that the LO predictions for A^(pt) are practically independent of (5 = ms/mc at pt >1 
GeV. 



C. Radiative Corrections 

In Ref. [pl| , radiative corrections to the SSA in heavy quark production by linearly polarized photons have 
been investigated in the soft-gluon approximation. The calculations show that the azimuthal asymmetry in both 
charm and bottom production is practically insensitive to the soft-gluon corrections at energies of the fixed target 
experiments. This implies that large soft-gluon contributions to the spin-dependent and unpolarized cross sections 



cancel each other in Eq. (1.3) with a good accuracy. One can assume that the same situation takes place for the 
azimuthal asymmetry in the decay lepton distribution. 

Our calculations of the factorizable NLO corrections to Afj,{pT) show that this is really the case. We have computed 



both spin-dependent and unpolarized differential distributions (2^) of the heavy-quark photoproduction at NLO 
to the next-to-leading logarithmic accuracy. The NLO corrections to the width of the heavy quark semileptonic 
decays are known exactly [p^|j20|]. We have found that radiative corrections to the leptonic SSA, ^^(pt), in the 



reaction (L5) are of the order of 1% in the SLAG kinematics. 

Two main reasons are responsible for perturbativc stability of the quantity yl^(pT'). First, radiative corrections 
to the SSA in heavy quark production are small ||ll[]. Second, the ratio I^^'^ {x) / 1^°'^'^ (x) is a constant practically 
at all X, except for a narrow endpoint region a; w 1 (Note that /^°''"(x) is the LO invariant width of the 

semileptonic decay c I'^viXq given by ( ^.8[ ) while I^^'^{x) is the corresponding NLO one.) The details of our 
NLO analysis are too long to be presented here and will be given in a separate publication |15|. 



III. NONPERTURBATIVE CONTRIBUTIONS 



Let us discuss how the pQGD predictions for single spin asymmetry are affected by nonperturbative contributions 
due to the intrinsic transverse motion of the gluon and the hadronization of the produced heavy quark. Because of 
the relatively low c-quark mass, these contributions are especially important in the description of the cross section 
for charmed particle production At the same time, our analysis shows that nonperturbative corrections to the 
single spin asymmetry are not large. 



A. Fragmentation 

Hadronization effects in heavy flavor production arc usually modeled with the help of the Peterson fragmentation 
fmiction pH , 

D{v) 2 , (3.1) 

where is a normalization factor and Sd — 0.06 in the case of a ZJ-meson production. The double differential 
distribution of a heavy meson has the form 

iPur)=[^Diy),,^ipnT/y), (3.2) 



dpDTdipQ J y dpQTdLpQ 
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where dpQ3°?yQ inclusive spectrum of the reaction (LI) and ym — 



Our calculations of the asymmetry in a D-meson production at LO with and without the Peterson fragmentation 
effect are presented in Fig. ^ by dotted and solid curves, respectively. For pdt > 1 GeV the fragmentation 
corrections to Ac{pt) are less than 10%. 

Analogous corrections to the asymmetry in the decay lepton azimuthal distribution, A^(j)t), are given in Fig. 



One can see that the effect of the fragmentation function (3.1) is practically negligible in the whole region of piT 



B. \lt Smearing 

To introduce fc^ degrees of freedom, kg ~ zk^ + kx^ one extends the integral over the parton distribution function 



in Eq. ( p. 151) to fcr-space, 

dzg{z, jip) dzd'^kTf (kr^ g{z, fip). (3-3) 
The transverse momentum distribution, / (kr^ , is usually taken to be a Gaussian: 



c 



In practice, an analytic treatment of kx effects is usually used. According to p2|, the fcx-smeared differential cross 



section of the process (1.1) is a 2-dimensional convolution: 



dpQTd(pQ ^^'^'^^ J d kr ^^^2^ dpgrdipg y'^^ 2^^ ' ' ^^'^"^ 



The factor i in front of kr in the r.h.s. of Eq. ( |3.5| ) reflects the fact that the heavy quark carries away about one 
half of the initial energy in the reaction (^]|) . 

Values of the fc^-kick corrections to the asymmetry in the charm production, Ac{pt), are shown in Fig. |^ by 
dashed ((fef^) — 0.5 GeV^) and dash-dotted {{k^} = 1 GeV^) curves. One can sec that /cT-smearing is important 
only in the region of relatively low pqt < rric- Note also that the fragmentation and A;T-kick effects practically 
cancel each other in the case of (fc|n) = 0.5 GeV^. 

Corresponding calculations for the case of the lepton asymmetry are presented in Fig. |^. It is seen that Afj_{pT) 
is affected by fcr-corrections systematically less than Ac{pt)- 

C. Fermi Motion 

The third type of nonperturbative corrections we considered can be associated with the motion of the heavy 
quark inside the produced (and decaying) hadron; they are commonly referred to as Fermi motion. These effects 
are included in the heavy-quark expansion by resumming an infinite set of leading-twist corrections into a shape 
function F (fc+), where k^ = ko + k^ is the positive light cone component of the heavy quark residual momentum 
inside the heavy meson [^|4j . The physical decay distributions, dFg^^'' , are obtained from a convolution of parton 
model spectra, drgj^\ with this function. This convolution is such that, in the perturbative formula for the decay 
distribution, the heavy quark mass is replaced by the momentum-dependent mass m* = mQ+ fc_|_ [p5| : 

rl3r(H) p d^r'*^^ / M \ 

E,-r^{MH,XH)= / '3 {m*,XH^)F{m*)dm*. (3.6) 
d-^P£ J d-'pi \ m* J 



xhMh 



The parton-level width, d'^r^l^'' (mg, x), in the r.h.s of Eq.(3.6) is defined by Eq. (|2.8|), Mh is the heavy meson mass 



and the hadron level variable xh is related to a; = 2[p£ ■pQ)/mQ by xh — xtuq/Mh- 
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Some of the properties of F {k+) are known. The first three moments satisfy Aq = 1, Ai = and A2 = ^mJj where 
/ij is the average momentum squared of the heavy quark inside the heavy meson and the moments are defined by 

A„ = (fc!^> = /_„^dfc+fc!^F(A;+) with K = Mh - mq. In our calculations we use a simple two-parametric model 

11: 

F(fc+)-iVe(i+'')''(l-r?r; ry=^, (3.7) 

where the condition Aq = 1 fixes the normalization factor N and the parameter a is related to the second moment 
as A2 = i/i^ = A^/(l + a). _ 

Our predictions for the asymmetry A^(pt) at different values of and A are give in Fig. ^ One can see that 
the bound state effects due to the Fermi motion of the heavy quark inside the heavy hadron are small in the whole 
region of pix- Taking into account that A^{pt) is also practically independent of rus at piT > 1 GeV, we conclude 
that our predictions for the SSA in secondary lepton distribution arc insensitive to the theoretical uncertainties in 
description of the charm semileptonic decays 

IV. EXPERIMENTAL CONSIDERATIONS 

The conditions needed to measure the SSA in open charm photoproduction are: 

• a high intensity photon beam with energy above 30 GeV, 

• a high degree of linear polarization, and 

• a large acceptance spectrometer. 

These conditions are similar to those of the approved E161 experiment |Q at SLAG, in which circularly polarized 
photons impinging on longitudinally polarized nucleons will be used to study AG{x) via open charm production. 
To measure the SSA under consideration, target polarization is not needed, so the microwaves that polarize the 
planned LiD target can be turned off. Linearly polarized photons can be obtained with an appropriate choice of the 
diamond used to produce coherent bremsstrahlung. Two orientations of linear polarization are obtained by rotating 
the diamond. The SSA can be measured using high pT muons in the E161 spectrometer, with cuts that minimize 
backgrounds. 

A. Photon Beam 

The most important difference from E161 is the change from circular to linear polarization. The E161 experiment 
plans to use coherent bremsstrahlung from a diamond target to produce polarized photons. Longitudinally polarized 
electrons are needed for circular polarization, which is largest for large values of E^/E,,, where E^ is the incident 
electron beam energy, and Ery is the resulting photon beam energy. To obtain linear polarization, the electrons can 
be unpolarized, and the extent of linear polarization is largest for low values of E^/E^. For a fixed photon energy, 
it is therefore essential to use the highest possible value of E^, which is about 50 GeV at SLAG. Traditionally, 
thin diamonds and a high degree of photon coUimation are used to enhance the ratio of coherent to incoherent 
photons, and hence the effective polarization. To obtain a sufficient flux of photons at SLAG, it is necessary to use 
a relatively thick diamond (we choose 1.5 mm), and a modest coUimation angle of 35 fii, which is several times the 
characteristic bremsstrahlung angle m^/ E^ « 10 /zr. The photon intensity and the degree of linear polarization for 
a diamond orientation that produces a primary peak from the (022) inverse lattice point at 35 GeV are shown in 
Fig. ^ for the emittance and maximum intensity conditions of the SLAG electron beam. Galculations were done 
using a Monte Garlo simulation based on the formulas of Ref. |^ . The flux in the region 30 to 35 GeV is about 
10^" photons/sec. The linear polarization peaks at about 0.4, with an average value of about 0.2 in the 30 to 35 
GeV region. Gontributions from the higher energy peaks (from the (044), (066), etc. inverse lattice points) will be 
small because both the flux and polarization are considerably reduced compared to the primary coherent peak. 
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B. Open Charm Detection and Backgrounds 



Because of the poor duty factor at SLAC, simulations made for E161 have shown that the best signal-to- 
background ratio is obtained by tagging open charm with one or more high transverse momentum muons. In 
the case of the SSA, the asymmetry rises rapidly with px, then becomes roughly constant for pT > I GeV. All 
background processes except J ftp decay decrease more rapidly with pT than open charm (as discussed in more detail 
below), therefore measuring at the highest possible pt will optimize signals over backgrounds. However, the open 
charm cross section decreases with px, which reduces statistical accuracy. The best compromise is for px values 
centered at 1 GeV. 

The most significant background process generating single muon with high px is the decay of pions and kaons. 
This was simulated using PYTHIA for typical photoproduction events that were allowed to decay in a block 
of copper placed close to the target. Studies using GEANT indicate that an effective path length before pions and 
kaons are absorbed is 50 cm of copper, so this was used as the effective decay length. The rates from tt, K decays 
are shown as the diamonds in Fig. |l^ for any sign muon with 7.5 < £"* < 10 GeV. They are compared with the rate 
of muons predicted from open charm decay as calculated using the photon-gluon fusion process in PYTHIA, scaled 
by a factor of two to match experimental data. The tt, K decay background is lower for negative muons than for 
positive muons by typically 50%, so one sign will be easier to measure than the other. The rate and SSA for tt, K 
decay muons can be measured by moving the copper absorber a significant distance from the target, thus doubling 
or tripling the decay rates. Measurements can also be made for 5 < -E* < 7.5 GeV, with somewhat larger relative 
backgrounds from Tr,K decays. In the interesting range 0.8 < px < 1.2 GeV, tt and K decays become dominant 
below E* w 4 GeV, making experimental measurements difficult in this energy region. 

Another significant process is wide-angle muon pair production from the target nuclei (Li and D). The rates 
calculated using Ref. ||2^ are shown as the crosses in Fig. |l^. The rates are for events in which only one of the two 
muons from the Bethe-Heitler pair is detected in the proposed E161 spectrometer. Events in which both muons 
would be detected can be vetoed in the data analysis, reducing the Bethe-Heitler background by about a factor 
of two. The remaining background is well below the open charm rate. It can be calculated quite accurately, and 
reliable corrections applied to the data. 

Backgrounds from the decay of vector mesons are relatively unimportant for px ~ 1 GeV, and can be measured 
using the approximately 50% of events in which both muons are detected in the spectrometer. 

The final background is from J/iIj decay, for which the typical px peaks around 1.5 GeV. Fortunately, most of 
the time the second muon is detected in the E161 spectrometer, so the background only becomes large for px > 1.3 
GeV and E'^ > 10 GeV, and can be measured using the large sample of events in which both muons are detected. 

In summary, in the optimal case of negative muons with 7.5 < < 10 GeV and px > 0.9 GeV, the signal-to- 
background ratio is estimated to be about 3:1, leading to an effective dilution factor of about 1.3, equivalent to a 
factor of 1.6 more running time than is needed in the absence of backgrounds. This estimate relies crucially on 
the open charm cross section as a function of pr, which is poorly known at present (factor of two uncertainties). 
Thus, the relative background situation could be better or worse, depending on what open charm cross sections are 
measured in the experiment. 

C. Projected Errors 

The cross section for producing negative muons from open charm decays with 7.5 < E* < 10 GeV and px > 0.9 
GeV is approximately 0.2 to 0.5 nb at E.y — 35 GeV, as estimated using the PYTHIA Monte Carlo. If we assume 
the E161 luminosity, i.e. a flux of 10^° photons per second, a 1.6 gm/cm^ target, and a spectrometer acceptance 
of about 50% for the events of interest, the open charm rate is of order 100,000 events per day. For an average 
linear polarization of 0.2, a statistical error of about 0.02 on the SSA would be obtained in one day of running. 
This error bar is sufficient to make a good test of the QCD predictions, and further improvements to the statistical 
error would be soon be lost in the estimated 10% (relative) systematics error, dominated by the uncertainty in the 
linear polarization. With about a factor of four larger running time, measurements with an error of about 0.02 
could be separately made of positive and negative muons, and in several bins of i?* spanning 5 < E* < 20 GeV. 
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This would be very useful in testing the assumed dominance of the photon-gluon fusion mechanism. A difference 
in positive and negative muon results might indicate contributions from soft associated production processes, while 
high-Xi? behavior of the SSA can be sensitive to diffractive contributions |]lO|. Approximately the same running 
time would be needed for 25 GeV incident photons as for Ej = 35 GeV, because the lower charm cross section and 
higher backgrounds would be compensated for by higher photon linear polarization. 



V. CONCLUSION 



In this paper we analyze the possibility to measure the SSA in open charm photoproduction in the E160/E161 
experiments at SLAG where a coherent bremsstrahlung beam of linearly polarized photons with energies up to 40 
GeV will be available soon. In these experiments, charm production will be investigated with the help of inclusive 
spectra of secondary muons. The SSA transferred from the decaying c-quark to the decay muon is predicted to 
be large for SLAG kinematics; the ratio Ag{pT) / Ac{pt) is about 90% at pr >1 GeV. Our calculations show that 
the SSA in decay lepton distribution preserves all remarkable properties of the SSA in heavy flavor production: 
it is stable, both perturbatively and parametrically, and practically insensitive to nonperturbative contributions 
due to the gluon transverse motion in the target and heavy quark fragmentation. We have also found that QCD 
predictions for A£{pt) depend weekly on theoretical uncertainties in the charm semileptonic decays. We conclude 
that measurements of Ai{pt) in the E160/E161 experiments would provide a good test of pQCD applicability to 
open charm production. 

Note also that because of the low c-quark mass, the power corrections to the charm production can be essential. 
As was shown in |l^], the pQGD and Regge approaches lead to strongly different predictions for the single spin 
asymmetry in the region of low px and large Feynman xp- Data on the pt- and Xi^-distributions of the SSA in 
_D-meson photoproduction could make it possible to discriminate between these mechanisms. 
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FIG. 1. Comparison of the QCD LO pedictions for ^^(pt) and Ac{pt)- 
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FIG. 2. Dependence of the leptonic SSA, A^(pt), on the cutoff parameter cut- 



13 




0.5 1 1.5 2 2.5 

Pt (GeV) 



FIG. 3. SSA, A^{pt), in different intervals of the lepton energy E*. 
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FIG. 4. Dependence of the SSA in chaxm production, Ac{pt), on the c-quark mass, rric. 
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FIG. 5. Dependence of the leptonic SSA, Afi{pT), on the charm and strange quark mass, rric and m, 
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FIG. 6. SSA in a D-meson production; the QCD LO predictions with and without the inclusion of the kr smearing and 
Peterson fragmentation effects. 
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FIG. 7. SSA, Ay,{pT), in the decay lepton distribution; the QCD LO predictions with and without the inclusion of the kr 
smearing and Peterson fragmentation effects. 
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FIG. 8. Fermi motion corrections to the leptonic SSA, Ay,(pT), at different values of //^ and A. 
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FIG. 9. The upper panel shows the predicted relative intensity of photons (intensity is energy-weighted flux) for coherent 
bremsstrahlung with a primary coherent peak at E-y = 35 GeV. The lower panel shows the calculated linear polarization 
under conditions typical of SLAG E161. 
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FIG. 10. Relative number of counts from open charm (circles), pion and kaon decays (diamonds), nuclear pair production 
(crosses), light vector meson decays (stars), and J ftp decays (boxes), for 7.5 < i?* < 10 GeV as a function of pt for 
30 < E-f < 35 GeV. The count rates correspond to about one minute under E161 conditions. 
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